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^ ; Abstract 

Cylinder-shaped perfect lens deduced from the coordinate transformation method is proposed. 
<D ■ 

The previously reported perfect slab lens is noticed to be a limiting form of the cylindrical lens 
when the inner radius approaches infinity with respect to the lens thickness. Connaturality between 
a cylindrical lens and a slab lens is affirmed by comparing their eigenfield transfer functions. We 
numerically confirm the subwavelength focusing capability of such a cylindrical lens with consider- 
ation of material imperfection. Compared to a slab lens, a cylindrical lens has several advantages, 

O ' including finiteness in cross-section, and ability in lensing with magnification or demagnification. 
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£>y Immediate applications of such a cylindrical lens can be in high-resolution imaging and lithography 

Qh technologies. In addition, its invisibility property suggests that it may be valuable for non-invasive 

—i. \ electromagnetic probing. 
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The recent works on invisibility cloaking devices [l|, |2| have triggered a widespread interest 
on design of functioning electromagnetic (EM) devices based on the coordinate transforma- 
tion approach. In an effort to unify a range of EM "meta-phenomena" , Leonhardt and 
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Philbin point out that a perfect slab lens made of negative index material (NIM) [3j can be 
interpreted as a result of a coordinate transformation that maps a single region in virtual 
EM space to multiple regions in physical space [4j . A recent extension of this concept for de- 
signing slab lenses with additional operation functions, such as image translation, rotation, 
and magnification etc, has been presented in Ref. [5j. In this paper we report deployment 
of a spatial mapping in cylindrical coordinate for designing an annular structure which can 
be considered as a cylindrical analog of the previously much debated perfect slab lens. Such 
a cylindrical lens overcomes two shortcomings of a slab lens. First, its body has a finite 
cross-sectional size. No structural truncation in cross-section is necessary for its physical 
implementation [6|]. Second, it can form magnified or demagnified image. These advan- 
tages can make such a lens more favourable across many scenarios. It should be noticed 
that Pendry has previously proposed a version of perfect cylindrical lens which is based on 
rolling up a perfect slab lens existing in the virtual EM space 2|- However, some material 
parameters of Pendry's cylindrical lens are gradient over the entire domain; in addition, the 
gradient parameters approach infinity in the interior of the lens. Hence, a direct forward 
theoretical confirmation of that lens is seemingly impossible, not to mention the difficulty 
in its realization. In contrast, the cylindrical lens under current study has all finite material 
parameters, and an appropriate design can leave the exterior region unaltered as free space. 
Therefore, the final device is much more favourable from both realization and application 
points of view. 

The proposed cylindrical perfect lens is realized by taking a coordinate transformation 
from a virtual free EM space (r', 9', z') to the physical space (r, 9, z), both denoted in 
cylindrical coordinate. To achieve imaging in radial direction, we map the coordinates as 
r' — /(r), 9' — 9 and z' = z. One simple example of /(r), as illustrated in Fig. [U(a), takes 
the form of 

sa+(l— s)b -r . 

* — r, it r < a; 
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sr + (l-s)b, if a < r < b; (1) 

r, if r > b, 



where a and b are the inner and outer radii of the deduced cylindrical lens (hence thickness 
d = b — a) , and s is the negative slope deployed in the mapping function. The continuity 
of the radial mapping function is required in order to achieve impedance-matched material 
interfaces. The negative slope at a < r < b is responsible for transforming a single region in 
virtual space to multiple regions in physical space. While the slope can be of any negative 
value, in this paper we focus mainly on a slope of s = — 1. In Eq. [TJ the space outside 
the lens body is kept unaltered, which allows the lens structure to work in a free space 
environment. 

Coordinate transformation projects field in EM space directly to physical space in ac- 
cordance to the form-invariance property of the Maxwell equations. Consider a line current 
source which is positioned at r' = ^^ and oriented along z' axis in virtual space. Its 
electric radiation field E' z , (a zeroth-order Hankel function) can be mapped to the physical 
space via the relation E z {r) = E' z ,{r'). The resulting field distribution in physical space, 
illustrated in Fig. Hfb) by its real part, shows as if there are three current sources located 
at r = ^p (outside), r = ^ (within), and r = 3 ^2 a (inside). The equivalence of such 
multiple points in physical space implies that the device is able to duplicate perfect copies 
of a source [J, [8(. The lens can image either in inward mode (object is placed outside and its 
image is formed inside) or reversely in outward mode with comparable performances. Given 
the mapping function described by Eq. [Q the object needs to be placed within a certain 
distance to the lens' surfaces: outside the lens, the effective region is b < r < {2b — a); and 
inside -^r— < r < a. 

26— a — u 

The mapping function in Eq. [J can be readily used to understand the so-called hyper- 



lenses proposed in Refs. [9|, llC], llll] under the context of cylindrical structure. These designs 
were proposed to solve the impedance mismatch problem of the multilayered cylindrical 
lenses theoretically analyzed in [12|, ll3( and demonstrated in [14J, ll5[ . Under the framework 
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of transformation optics, the designs in [9|, llQ|, LL3J are inherently "less per fect" than the cylin- 
drical lens proposed in this study. The design method proposed in jljj can be interpreted 
as using an identical mapping function as in Eq. dj but with a zero slope for a < r < b. 
For such a lens, perfect imaging happens only when the object is placed exactly at r = a, 
and a corresponding image at r = b. Notice also that a coordinate transformation with a 
zero-sloped mapping function leads to unphysical extreme material parameters. While in 
[9J, the design of a hyperlens can again be interpreted as using Eq. [TJbut with a close-to- zero 




FIG. 1: Schematic illustration of the principle for idealized perfect cylindrical lensing. (a) Spatial 
mapping function, (b) Imaging of a line current source by a perfect cylindrical lens. Colormap 
shows the real part of the E z field. 

positive slope for a < r < b. No perfect imaging can be achieved in this scenario. The lens 
relays an object's near field at r = a to r = b with minimal decaying, owing to the adjacency 
of r = a and r = b surfaces in virtual space. 

According to the coordinate transformation described by Eq. [U we derive material pa- 
rameters for the designed cylindrical lens body (a < r < b) as 

r r — 2b 



r -2b 

e r — \l r — , €g — fig 



r -2b' 



e 2 = H z 



(2) 



and for interior region r < a as 

€ r = fi r = l,£g = fig = l,e z = fl z 



2b -a 



(3) 



The lens has negative anisotropic permittivities and permeabilities. At r = b, the lens has 
all material parameters valued at —1, which is perfectly matched to air; at r = a, it has 



e r = fi r = s -^-, eg = fig = -rzohi an d e z = l^z — ^"jr^ which matches to the interior material. 
EM field in the interior region is compressed in wavelength as compared to in free space 
due to its larger e z and fi z values. This contributes to image demagnification for inward 
lensing operation, or magnification vice versa. It is worth mentioning that when the lens 
thickness d is very small compared to a, the lens' material parameters at r = a approach to 
— 1, whereas e z and fi z values of the interior material approach to 1. That is, the cylindrical 
lens physically becomes a previously reported prefect slab lens [3|. 

The perfect cylindrical lensing phenomenon predicted by coordinate transformation, as 
illustrated in Fig. [Qb), however should be verified using full- wave analyses. One factor for 
the necessity of such explicit verification is that coordinate transformation has presumed 
the time-harmonic steady state. In a previous time-domain analysis for a slab lens, it was 
noticed that re-construction of a Fourier field component at the other side of the lens can 



take unpractically long time when the slab index approaches to ideal value, i.e. -1 16]. 



This suggests that perfect imaging by a slab lens is rather a "wishful" event which relies on 
infinite happening time. Nevertheless, subwavelength imaging can be realized by slab lenses 
with imperfect materials [16] , Il7| ■ Our frequency- domain finite element calculation showed 
no convergence for a cylindrical lens with ideal parameters as specified by Eq. |2j which 
reflects the same singular problem of the lens. In the following analyses, we examine the 
performance of cylindrical lenses with the presence of material perturbations. 

For ease of clarification, we restrict our studies to in-plane propagation case (wave vector 
in r9 plane) with electric field solely polarized in z direction (E z polarization). It should 
however be kept in mind that coordinate transformation promises focusing in a full 3D 
scenario with an arbitrary wave polarization. The relevant material parameters hence are 
e z , /i r , and fig. We append a common factor g, with g = 1 + 5, to each of the lens' 
ideal material parameters. S is referred to as material perturbation, and its imaginary part 
represents material loss. Material perturbation, especially loss, is inevitable for resonance- 
based negative permittivity or permeability metamaterials. For convenience, we denote 
layers in the three- layered system from inside to outside as layer 1, 2 and 3 respectively. 



General field solutions in three layers, in terms of E z , are 

E l = J2^-' y - J ' 
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in 



(4) 
(5) 
(6) 



where k$ is free-space wave number, m is the azimuthal momentum number, and srf % and 
SS l {i = 1,2,3) are coefficients. We take exp(— jut) time harmonic dependence. Therefore 
the Hankel function of the first kind Hm in Eqs. HH6] represents an outgoing cylindrical 
wave. Unlike in a planar system where purely evanescent or propagating eigenwaves exist, 
all eigenwaves in the cylindrical lens system are propagating. Despite the dissimilarity in 
eigenwaves, we show in the following that superlensing mechanisms of two types of lens can 
be bridged. 
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FIG. 2: Comparison of transfer functions for a cylindrical lens and a slab lens of the same thickness. 
Three different object positions are shown: for cylindrical lens r 0j i = ^§rEi~ ? r o,2 = 3 ^,_ 2° ? r o,3 = 
86-la ' anc ^ f° r s ^ a b l ens x o,i = — X' x °> 2 = — ^' x °>3 = — T - Two sets of object positions are 
equivalent when o and b are much larger than d. For the cylindrical lens the corresponding angular 
momentum number is marked on the secondary x axis. 



First, we consider a cylindrical lens whose inner radius a is much greater than both 
its thickness d and the operating wavelength. The system resembles a slab lens in both 



geometry and material parameters. Our objective is to show how the imaging performances 
of such a cylindrical lens and a slab lens of the same thickness can be quantified in a 
unified manner. For a slab lens, the transfer function, which describes the ratio of field 
strengths at image and source positions for each plane wave component, has been used for 



describing its imaging capability [16|, lljj, ll8| . Consider a slab lens with thickness d = A/2 
which is placed in yz plane and centered at x = 0. Likewise, we assume that waves are 
propagating in xy plane with pure E z polarization. The slab lens has material parameters 
of e = fi = —g = —1 — 5, with 5 = 0.0001 — O.OOOOli. Its transfer functions, as characterized 
by normalized k y are shown in Fig. [2] (solid curves) for three different object positions. For 
this particular material perturbation introduced, restoration of evanescent waves ik y > k ) 
at the image positions (always 2d apart from the source) has been realized up to a certain 
k y . The transfer functions are characterized by a resonance before the transmission damps 
heavily. The free-space transfer function from the source position to image position (without 
the slab lens) is also included in Fig. [2l The evanescent waves can barely reach the image 
position. For a cylindrical lens, one may first be puzzled by the absence of decomposable 
wavevectors in the general solutions (Eqs. HUE]). In fact, the azimuthal momentum number m 
determines the pace of field variation along the surface-tangent (8) direction. The azimuthal 
momentum number can be converted to azimuthal wave number as kg = m/ri, where r, 
is the radial image position [19j. In order to obtain sharp image especially in 9 direction, 
existence of high-order cylindrical wave components at the image position is essential. Hence, 
transmission of cylindrical waves through a cylindrical lens provides vital information for 
the lens' imaging capability. Here we consider a cylindrical lens with d = A/2 and a = 10A. 
5 is kept the same as in the slab case. The cylindrical wave transfer functions for outward 
lensing operation, plotted in terms of normalized kg, are superimposed in Fig. [21 Three 
different object positions have been studied in Fig. [2j The transfer function records the 
ratio of field strengths at the image and source positions at different cylindrical wave order 
number m. For each m, the field ratio is derived through matching of the tangential fields 
at both r = a and r = b interfaces. Referring to Eqs. \ME\, for outward imaging operation, 
the task is to find the coefficient &P from a known value of SB 1 ; while for inward imaging 
operation, the task is to find the coefficient stf x from a known value of &/ 3 . Unlike for a 
slab lens case, the transfer function for a cylindrical lens is too lengthy to be expressed here 
owing to the presence of transcendental Bessel functions. For consistency, the free-space 



cylindrical wave transfer function is also shown. From Fig. [21 we notice excellent agreement 
between the transfer functions for cylindrical and slab lenses. With the same thickness 
and material perturbation, the two types of lenses are therefore equivalent in their imaging 
capabilities. Both of them can restore a spectrum of waves {k{e, y } > 1) that were lost in 
the absence of the lenses. Figure H] provides independent verification of the sub- wavelength 
imaging capability of a cylindrical lens. 
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FIG. 3: Effect of a on cylindrical wave transfer function. Inset gives the zoom- in view of the 
transfer function with a = 1.5A (inward). 

Next, we look into the imaging performance of a cylindrical lens as its inner radius a 
decreases. At a smaller a value, the physical problem becomes more distinct as compared to 
a slab lens. It is worth noting that, when the object and image points under examination are 
very close to the origin of cylindrical system, a cylindrical wave can appear to be a strong 
mixture of both propagating and evanescent waves. This can be observed by comparing 
the free-space transfer functions in slab and cylindrical systems shown in Fig. [3j Having 
known this, one would naturally expect to see the distinct character of a cylindrical lens 
as reflected in its transfer function. The result however turns out to be a bit surprising. 
The cylindrical wave transfer function is hardly affected when its inner radius is reduced 
even close to wavelength. In our case study, we fix thickness d = 0.25A and decrease a 
from infinity to 0.5A. Material perturbation is again held at 8 = 0.0001 — O.OOOOli. Besides 
showing the effect of reducing a on the transfer function, we also take this opportunity to 
examine the difference in performances between two operation modes of the same lens, i.e. 
inward imaging and outward imaging. For each operation mode, we start from the slab 
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lens case (a = oo), and decrease a to 1.5A and then to 0.5A. The variation of the transfer 
function is shown in Fig. [31 Notice that for inward operation, kg is normalized with respect 
to the wavenumber in layer 1, i.e. k\ = ko 2 ^ 11 ^. Object position for inward operation is at 
(36 — o)/2, and for outward operation (3ab — a 2 ) /(4b — 2a). Figure [3] shows that the transfer 
function shifts only very slightly when a decreases from infinity to 0.5A. Therefore it can 
be concluded that, with the same thickness and material perturbation, a cylindrical lens, 
even at an overall scale comparable to wavelength, is as effective as a slab lens for achieving 
subwavelength imaging. This conclusion is valid for either inward or outward operation 
mode. The sharp variations of the inward transfer functions at small ke values (see the inset 
in Fig. [3]) are due to near zero field values at both object and image positions. They are 
not caused by resonances. 

The inevitable material imperfectness imposes a rather stringent limitation on the lens 
shell thickness. For a slab lens, it has been noticed that a smaller thickness drastically 
improves the imaging performance (or equivalently, relaxes the requirement for material 



perfectness in order to achieve the same subwavelength imaging) 17j |. We have noticed the 
same for cylindrical lenses. This information can be obtained by comparing Fig. [2] and Fig. 
|3l In general, with a reasonable material perturbation, the thickness of a cylindrical lens 
should be kept less than a wavelength in order to achieve subwavelength resolution. Here 
we explicitly look into the effect of material perturbation on a relatively thin cylindrical 
superlens' performance. The cylindrical lens to be studied has fixed geometrical parameters 
as a = 0.5A and d = 0.25A. First the effect of material loss is examined by keeping the 
real part of 5, or 3ft(<5), at 0.0001, and increasing the imaginary part S(5) from 0.00001 up 
to 0.001. The variation of the transfer function is shown in Fig. |U From the figure, it is 
noticed that if Q(S) < 9?(<5), only the "finesse" of the resonance peak is affected. Higher 
loss suppresses the resonance. Only when S(<5) > 9?(<5), the transmission starts to drop 
at a significantly lower kg value, which therefore affects the image resolution. A material 
perturbation dominated by either its real part or imaginary part (but with around the same 
value) will result in similar transfer functions, with the exception that no resonance peak is 
present for the transfer function resulted from the latter perturbation. This is confirmed by 
two curves (5 = 0.0001 — O.OOli and 5 = 0.001— 0.0001i) in Fig. HI The imaging enhancement 
factor, defined as the maximum kg/ki at which the transmission from the object to image 
is close to unity is about 4.5 for 5 = 0.0001 - 0.001« and 6 for 5 = 0.0001 - O.OOOOli. The 



results suggest that cylindrical lenses under study are able to image with subwavelength 
resolution. When the perturbation is reduced further, the resonance peak will shift to an 
even larger kg value. In the limit, for an ideal lens the resonance happens at infinite kg. 
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FIG. 4: Effect of 8 on cylindrical wave transfer function. Azimuthal momentum number is marked 
on the secondary x axis. 



In Fig. [51 we show the radial \E Z \ field distribution of a high-order cylindrical wave 
(m = 15) for the lens studied in Fig. HI Particularly, we increase the loss perturbation and 
observe its effect on field profile. di(S) is kept at 0.0001, and 3(5) is changed from 0.00001 
to 0.005, 0.05, and 0.1. The study is for inward imaging scenario. The cylindrical wave, at 
the vicinity of the object or image position, resembles a purely evanescent wave. The field 
only becomes oscillating at a distance further away from the origin. With a very small loss, 
the incoming field from outside penetrates freely (through the formation of a so-called anti 
surface mode at the outer interface) into the lens body and excites a surface mode at the 
inner lens boundary. The field amplitude is therefore amplified. With a further decay in 
lens interior, the field amplitude at the image position is almost exactly the same as at the 
object position. Restoration of the cylindrical wave is realized. When the loss increases, a 
surface mode starts to build up at the outer boundary of the lens, while the strength of the 
surface mode at the inner boundary is gradually muted. At 3(5) = 0.01, the loss becomes 
so heavy that the first surface mode also weakens. The results shown here are very similar 
to those noticed for a slab lens 20J]. 

To complement the above transfer function analyses, here we numerically demonstrate the 
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FIG. 5: Radial |E 2 | field distributions as material loss increases. The vertical solid black lines 
denote lens boundaries; dashed red lines denote object and image positions. 

imaging performance of cylindrical lenses using the finite element method (FEM). The issue 
of coupled surface-mode resonance, which gives rise to the peaks in the transfer functions, is 
also addressed here. The first lens to be examined has parameters of a = 0.5A and d = 0.25A, 
and S = 0.0001 — O.OOOOli. We position a z-oriented line current source at r = (36 — a)/2. 
An image is expected at r, = (3ab — a 2 )/(4b — 2a) in lens interior. The real part of the steady 
state \E Z \ field is captured in Fig. [6]^a). An image is clearly re-constructed inside the lens. 
However, the most noticeable difference as compared to the idealized field distribution in Fig. 
[T]is that strong localized fields are present at the two surfaces of the lens. The maximum field 
values on the surfaces are as high as 6500 positively and -5000 negatively (in arbitrary unit), 
which are about 10 times as large as the value at the object position. There are altogether 
70 nodal lines along azimuthal direction, indicating the resonance happens around m = 35, 
which agrees with its corresponding transfer function in Fig. El The imaginary part of 
the field (not shown) has the same number of nodal lines, but with different positions as 
compared to the real part. Hence the surface wave is not standing, but traveling along 6 
direction. An examination of the Poynting vector plot confirms giant energy flow around 
the lens' surfaces (not shown). The huge amount of energy that the lens can trap suggests 
that it may require a certain amount of time for the system to reach its steady state, and 
in turn to form the image. Another consequence of the existence of heavy surface modes is 
that it may disallow an object from being placed very near or very far to the lens surface, 
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which adds further restriction on the effective lensing and detecting areas. From Fig. HI we 
know that a relatively lossy material can help to suppress the resonance peak in transfer 
function, though at the expense of a lower achievable resolution. We therefore carried out 
another FEM analysis for the same structure but with 5 = 0.0001 — O.OOli. The real part 
of the E z field is shown in Fig. El^b). Compared to Fig. E](a), complete suppression of the 
surface mode is seen at surface positions away from the line source. Still, surface mode is 
present close to the line source, but its strength has been reduced considerably compared to 
the previous calculation. For both cases, despite the localized surface modes, field outside 
the lens experiences almost no disturbance by the lens structure; it appears to be radiated 
directly from the line source. Therefore the lenses are almost invisible to an observer. The 
invisibility is due to the fact that the particular coordinate transformation deployed assures 
the whole structure bounded by r < b is electromagnetically equivalent to vacuum. The 
perturbations introduced here are insignificant for affecting the invisibility property. 
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FIG. 6: Inward imaging of a line current source, (a) 5 = 0.0001 — O.OOOOlz. Overvalued field are 
shown in saturated color as dark red (positive) or dark blue (negative), and contour lines at the 
levels of ±1000, ±2000, ±3000, ±4000, ±5000 and 6000 are marked, (b) 5 = 0.0001 - O.OOli. 
Contour lines at the levels of ±1000, ±1500, 2000, 2500, 3000 and 3500 are marked. Yellow crosses 
are object and image positions. Domain size: 2.5A x 2.5A. 

With the same two lenses, we numerically demonstrate their inward imaging of two 
parallel line current sources. Such numerical experiments provide direct evidence for the 
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maximum achievable resolution. Note that here we calculate the distance between two line 
sources or their images by considering they are located on a cylindrical plane (concentric to 
the lens surfaces). Two z-oriented line sources are placed at r = (36— o)/2, and azimuthally 
separated by a distance denoted as ALg. Their images are recorded while ALg is decreased 
gradually. The recorded electric field intensity at object and image planes for the low- loss 
and lossy lenses are shown by the panels in the first and second rows in Fig. [71 respectively. 
For the low-loss lens, the intensity at the image plane shows very clearly the presence of 
two peaks when ALg = 0.22A, which corresponds well to the field at the object plane. As 
ALg decreases, the peaks become less evident, and they are overtaken by a central intensity 
peak when ALg = 0.16A. The central peak is caused by the extension of a strong surface 
mode. It is however arguable that from the image obtained for ALg = 0.16 A, one still can 
tell the information of two line sources. The overall image, though distorted, can possibly be 
corrected through signal processing techniques, especially if one has the preknowledge of the 
distortion. For the lossy lens, the information of the two current sources is almost completely 
lost at ALg = 0.16A. However, the two line sources becomes identifiable once ALg reaches 
0.18 A. The recorded field intensities show much less ripples due to the suppression of surface 
modes. In all panels in Fig. [7J the blue curves represent the field intensities at the object and 
image planes while the lens is absent. Drastic decay in amplitude is noticed for the images 
recorded without a lens. Hence the restoration of evanescent EM wave by a cylindrical lens 
is obvious. Comparing the low-loss lens to the lossy one, we notice that the latter suffers 
less disturbance from surface mode. Due to this factor, its maximum achievable resolutions 
for both lenses are almost comparable. 

Lastly, we make a couple of comments on the performance of a cylindrical lens when 
some parameter choices for designing the lens are different from above analyses. First, we 
comment on the effect of negative slope value used for the spatial mapping function (Eq. 
[T]). In practice, one desires to image an object that is far away from the lens, rather than 
maneuvering the lens very close to the object for performing an imaging. Theoretically this 
can be achieved by an increase (negatively) in the negative slope value for the mapping 
function. However, for a fixed lens thickness, an increase in the slope compresses the field in 
the lens body, or equivalently increases proportionally the optical lengths for all cylindrical 
wave components as they travel across the lens body. That unfortunately amplifies the 
adverse effect of material loss to the lens' performance. We noticed that in order to achieve 
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FIG. 7: Imaging of two parallel line current sources. Panels in the first row are for the lens with 
5 = 0.0001 — 0.00001i, and those in the second row are for the lens with 5 = 0.0001 — O.OOli. For 
panels in both rows, from left to right, ALg values are respectively 0.16A, 0.18A, 0.20A, and 0.22A. 
Notice that the image has been stretched by a factor of (26 — a) /a due to the higher permittivity 
inside the lens. Red lines: with lens; blue lines: without lens. Thick lines: intensity at image plane; 
thin lines: intensity at object plane. 

the same resolution, the thickness of a lens should be kept inversely proportional to the 
negative slope, assuming the same material perturbation. Due to the necessity of a decrease 
in lens thickness for maintaining resolution, the object-to-image distance in fact shortens 
when the negative slope increases. Another note is that an increase in the negative slope gives 
rise to gradient lens material with higher anisotropic ratio. Therefore such a way of increasing 
the object distance outside the lens is not encouraged. Our second comment is regarding the 
material perturbation introduced. So far in our previous theoretical studies we have chosen 
9ft(5) to be positive, while 3(<5) is kept negative. Negative $s(5) is natural for denoting loss 
incurred by passive materials. It is however foreseeable that in practical implementations 
9?(<5) can easily be negative. Through calculations we have noticed the following situations 
for choosing a negative di(5). When |§R(5)| is relatively larger than |S(5)|, change in the 
sign of di(S) from positive to negative will affect the characteristics of the transfer function 
only around the resonance peak position. Noticeably, the resonance peak will shift slightly 



14 



to a larger k$ position, regardless of inward or outward operation mode. For example, the 
resonance peak of the transfer function for the lens with 5 = 0.001 — O.OOOli in Fig. H]will 
shift from kg/ki = 4.70 to 5.05 if we use 5 = —0.001 — O.OOOli. Under the scenario that 
1 9ft(<5) | is relatively smaller than |3(5)|, no significant change in transfer function can be 
noticed when we flip the sign of 9?(£). 

In conclusion, based on the coordinate transformation technique we have designed a type 
of cylindrical superlenses. The transfer function analyses show that such a cylindrical lens, 
even with a wavelength-scale cross-section size, can achieve similar subwavelength imaging 
resolution as compared to the previously reported negative index slab lens. Subwavelength 
imaging is demonstrated numerically for cylindrical lenses with imperfect material param- 
eters. Like its slab lens counterpart, the cylindrical lens experiences similar limitation as 
imposed by inevitable material imperfectness. Most noticeably, at a practical loss level, 
a cylindrical lens should have a subwavelength thickness in order to realize subwavelength 
focusing. Therefore, these cylindrical lenses are bound to be "near-sighted", a conclusion 
previously drawn for slab lens [18(. Besides high- resolution imaging, the proposed cylin- 
drical lens can be applied for lithography applications. Consider a single E z polarization 
operation. A properly designed cylindrical lens can embed in its interior a dielectric material 
with perfect impedance matching. Image can then be written on the surface of that dielec- 
tric material. The invisibility property of the cylindrical device can also be deployed for 
non-invasive electromagnetic probing. Although not analyzed in this paper explicitly, the 
coordinate transformation predicts that the same radial transformation as described in Eq. 
[TJwill yield a spherical superlens. With the current fast-developing metamaterial technology, 
which has already been used for demonstrating negative refraction |2ll . |22J, subwavelength 
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231 ]. and even conceptual electromagnetic cloaking [24] (also devised in Ref 



251 ] for operating at optical wavelength), etc., we should be able to realize fabrication of the 



proposed cylindrical superlens in near future. 
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